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Abstract 



We evaluate the 0(a) QED corrections to the recoil electron energy spectrum in 
the process ui + e — > Vi + e (+7), where (+7) indicates the possible emission of a 
photon and I — e, fi or r. The soft and hard bremsstrahlung differential cross sections 
are computed for an arbitrary value of the photon energy threshold. We also study 
the O(a) QED corrections to the differential cross section with respect to the total 
combined energy of the recoil electron and a possible accompanying photon. Their 
difference from the corrections to the electron spectrum is investigated. We discuss 
the relevance and applicability of both radiative corrections, emphasizing their role 
in the analysis of precise solar neutrino electron scattering experiments. 
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1 Introduction 



The QED corrections to the scattering process v e + e — > v e + e were studied long ago 
by Lee and Sirlin using the effective four-fermion V— A Lagrangian JEJ, and shortly 
afterwards Ram || extended their calculations by including hard photon emission. A 
few years later, 't Hooft computed the lowest order prediction to this differential cross 
section in the Standard Model (SM) ||. Since then, the radiative corrections to this 
process, which plays a fundamental role in the study of electroweak interactions, have 
been investigated by many authors, focusing on various aspects of the problem [f§-[n[ • 
't Hooft 's early SM predictions for u-e scattering were used by Bahcall to exam- 
ine the total cross section, energy spectrum and angular distribution of recoil elec- 



trons resulting from the scattering with solar neutrinos |15|. More recently, Bahcall 



Kamionkowski and Sirlin performed a detailed investigation of the radiative correc- 
tions to these recoil electron spectra and total cross sections ||16|| . Their results show 
the importance of these corrections for the analysis of precise solar v—e scattering ex- 
periments, particularly of those measuring the higher energy neutrinos that originate 
from 8 B decay. 

In this paper we study the 0(a) QED corrections to neutrino electron scattering 
in the SM, with contributions involving either neutral currents (as in the + e — ► 
u h,t + e process) or a combination of neutral and charged currents (as in the u e + e — > 
v e + e process). In this analysis we make the approximation of neglecting terms 
of 0(q 2 /Ml), where q 2 is the squared four-momentum transfer and M w is the W 
boson mass. Within this approximation, which is excellent for present experiments 
(\q 2 /M^\ ~ 1 when the electron recoil energy ~ 6 x 10 3 TeV!), the SM radiative 
corrections to these processes can be naturally divided into two classes. The first, 
which we will call "QED" corrections, consist of the photonic radiative corrections 
that would occur if the theory were a local four-fermion Fermi theory rather than 
a gauge theory mediated by vector bosons; the second, which we will refer to as 
the "electroweak" (EW) corrections, will be the remainder. The EW corrections have 
been studied by several authors |5|-p|, plj||i6|1 and are not discussed in the present paper. 
The split-up of the QED corrections is sensible as they form a finite (both infrared 
and ultraviolet) and gauge-independent subset of diagrams. We refer the reader to 
ref. [[L7j for a detailed study of this separation. 

The QED radiative corrections are due to both loop diagrams (virtual corrections) 
and to the bremsstrahlung radiation (real photons) accompanying the scattering pro- 
cess. Of course, only this combination of virtual and real photon corrections is free 
from infrared divergences. To order a, the bremsstrahlung events correspond to the 
inelastic process v\ + e — » v% + e + 7 (I — e, fi or r). Experimentally, bremsstrahlung 
events in which photons are too soft to be detected are counted as contributions to 
the elastic scattering V\ + e — > v x + e. The cross section for these events should be 
therefore added to the theoretical prediction of the elastic cross section, thus removing 
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its infrared divergence. 

We will divide the bremsstrahlung events into "soft" (hereafter SB) and "hard" 
(hereafter HB), according to the energy of the photon being respectively lower or 
higher than some specified threshold e. We should warn the reader that the words 
"soft" and "hard" may be slightly deceiving. Indeed, if e is large (small), the SB 
(HB) cross section will also include events with relatively high (low) energy photons. 
While calculations of both soft and hard bremsstrahlung are often performed under 
the assumption that e is a very small parameter, much smaller than the mass of the 
electron or its final momentum, we will also discuss results for the case in which e 
is an arbitrary parameter constrained only by the kinematics of the process. Indeed, 
the HB cross section (contrary to the SB one) is by itself, at least in principle, a 
physically measurable quantity for any kinematically allowed value of this threshold. 
All calculations have been carried out without neglecting the electron mass. 

The paper is organized as follows. In sect. 2 we discuss the lowest order prediction 
for the neutrino electron differential cross section, together with its soft bremsstrahlung 
and one-loop QED corrections. The hard bremsstrahlung recoil electron spectrum is 
examined in sect. 3. In the same section, this contribution is added to the virtual and 
soft ones to derive the total QED corrections. In sect. 4 we evaluate the spectrum of 
the total combined energy of the recoil electron and a possible accompanying photon 
emitted in the scattering process. We summarize our main results in sect. 5, dis- 
cussing their applicability and emphasizing their role in the analysis of solar neutrino 
electron scattering experiments. 



2 Virtual and Soft Photon Corrections 

The SM prediction for the elastic neutrino electron differential cross section is, in 
lowest -order and neglecting terms of 0(q 2 /M^) H], 
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where m is the electron mass, 
sin 2 6 



constant 
and sin 6 W 



it 



1.16637(1) x 1(T 5 GeV~ 2 is the Fermi coupling 
±1/2 (upper sign for u e , lower sign for u^ T ), q R = sin 2 ^, 



g L = sin u w ± i/z (upper sign tor z/ e , lower sign lor f M , T J, g R 
is the sine of the weak mixing angle. In this elastic process E, the electron 
recoil energy, ranges from m to E max = [m 2 + (2u + m) 2 }/[2(2u + m)], z = (E — m)/u 
and v is the incident neutrino energy in the frame of reference in which the electron 
is initially at rest. We will refer to the L, R and LR parts of an expression to indicate 
its terms proportional to g 2 , g\ and g L g R , respectively. For example, the R part of 
[da/dE] (eq. ©) is (2mG 2 /7r)^ 2 (1 - zf . 

According to the definition discussed earlier, the one-loop QED corrections to 
neutrino electron scattering consist of the photonic vertex corrections (together with 
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the diagrams involving the field renormalization of the electrons) computed with the 
local four-fermion Fermi Lagrangian. These corrections give rise to the following 
expression for the differential cross section: 



where 
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A is a small photon mass introduced to regularize the infrared divergence and I 
\/E 2 — rri 
by 



is the three-momentum of the electron. The dilogarithm Li 2 (x) is defined 



Li 2 (x) 



t 



The L part of eq. (|j) (with g L = 1) is identical to the formula for the one-loop 
photonic corrections to the v e + e — > v e + e differential cross section computed long 
ago in the pioneering work of Lee and Sirlin using the effective four-fermion Fermi 
V— A Lagrangian. The analogous formula for the reaction involving an anti-neutrino 
V e (rather than a neutrino v e ) can be found in the same article and coincides[| with 
the R part of eq. @ (with g R = 1). This identity is simply due to the fact that 
the cross section for antineutrinos in the local V— A theory is the same as that for 
neutrinos calculated with a V+A coupling. On the contrary, the LR part of eq. (@) 
has clearly no analogue in the V±A theory, but can be derived very easily once the 
L and R parts are known. 



1 With the exception of a minor typographical error in their eq. 22, where the square bracket 
multiplying J ra d should be squared. We thank Alberto Sirlin for confirming this point. 
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The v e + e — > z/ e + e + 7 differential cross section with emission of a soft photon 
was computed in ref. Jl|, once again by using the effective four-fermion Fermi V— A 
Lagrangian. It can be identified with the L part (with g t — 1) of the soft photon cor- 
rections to the tree level result in eq. ([1]). The L, R and LR parts of these corrections 
(with g L = g R = 1) are however identical, because the whole soft bremsstrahlung cross 
section is proportional to its lowest-order elastic prediction. We can therefore write 
the soft photon emission cross section in the following factorized form: 
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(For x G 1R, L(x) = — Re[Li 2 (x)].) This result is valid under the assumption that e, 
the maximum soft photon energy, is much smaller than m or the final momentum of 
the electron. As we mentioned earlier, in the next section we will discuss numerical 
results for the case in which e is an arbitrary parameter. 

The reader will notice that the sum \Vi(E) + I~,(E, e)] does not depend on A, the 
infrared regulator. Indeed, the infrared divergence of the virtual corrections (eq. (0)) 
is canceled by that arising from the soft photon emission (eq. (H)). 



3 Hard Bremsstrahlung and Total QED Correc- 
tions to the Final Electron Spectrum 

The SM prediction for the differential neutrino electron cross section 

i/, + e->^ + e(+7), (8) 

where (+7) indicates the possible emission of a photon, can be cast, up to corrections 
of 0(a), in the following form: 
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(We remind the reader that terms of 0(q 2 /M^) are neglected throughout this paper.) 
The deviations of the functions g L (E) and g R (E) from the lowest-order values g L and 
g R reflect the effect of the electroweak corrections, which have been studied by several 



authors |8|, p~0| , p~6f| . (See ref. JT6| for simple numerical results.) 

The functions f x (E,u) (X = L,R or LR) describe the QED effects (real and 
virtual photons). For simplicity of notation their v dependence will be dropped in the 
following. Each of these functions is the sum of virtual (V), soft (SB) and hard (HB) 
corrections, 

f x (E) = f v x {E) + f x B (E, e) + f« B (E, e). (10) 

(We remind the reader that we have defined the bremsstrahlung events as "soft" or 
"hard" according to the energy of the photon being respectively smaller or higher 
than a specified threshold e.) The analytic expressions for f x (E) and f x B (E,e) can 
be immediately read from eqs. (§) and (g) respectively (the latter being valid only in 
the small e limit) and their sums, which are infrared-finite, will be denoted by 

K s (E,e) = fZ(E!) + fZ a (E,e). (11) 

Analytic expressions from which one can obtain f BB (E,e) and f" B (E,e) were cal- 
culated long ago by Ram ||. Although these results were obtained in the small e 
approximation, the formulae are nonetheless long and complicated and we will only 
plot the results for specific values of v and e. The function f BB (E,e) has not been 
previously calculated. 

We created BC, a combined Mathematica-FORTRAN codeQto compute the f x B (E, e) 
functions for arbitrary positive values of the parameter e up to the kinematic limit 
v (v, the incident neutrino energy in the laboratory system, is also the maximum 
possible energy of the emitted photon). We first computed the transition proba- 
bility for the bremsstrahlung process, averaged over the initial electron spins and 
summed over the polarizations of the final electron and photon. We then removed the 
energy-momentum conserving 8 function in the three-body phase-space integral and 
specialized the result to the laboratory frame of reference. For a given initial neutrino 
energy, the five independent variables describing the final state were chosen to be the 
four angular variables of the final electron and photon, plus the electron recoil energy. 
In order to compute the HB functions f x B (E, e) we finally imposed the condition that 
the photon energy u should be larger than the threshold e (without assuming e to be 
small). This constraint required a detailed analysis of the kinematically allowed ranges 
of variability of the chosen phase-space variables. The last integrations over the an- 
gular variables were then performed numerically using the Monte Carlo method [ 20f 



and demanding a 0.1% relative accuracy. (Reminding the reader that X indicates L, 
R or LR, we note that this uncertainty in the computation of a function f" B (E,e) 

2 The code BC, available upon request, computes all QED corrections discussed in this paper. It 
uses the Mathematica package FeynCalc fy| and the FORTRAN code VEGAS po|. 
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produces an extremely tiny relative error (a/ir)f^ B (E, e) x 0.1% in the corresponding 
X part of the differential cross section in eq. (|9|). This high level of accuracy was 
useful for internal numerical checks and for the comparisons, in the small e limit, with 
Ram's results.) All calculations have been carried out without neglecting the electron 
mass. 

In figs. 1 and 2 the functions f BB (E, e) are respectively plotted for v = 0.862 MeV 
and v = 10 MeV, setting the threshold e to several different values. These two values 
of v were chosen for their relevance in the study of solar neutrinos: v = 0.862 MeV is 
the energy of the monochromatic neutrinos produced by electron capture on 7 Be in 
the solar interior, while v = 10 MeV belongs to the continuous energy spectrum of 
the solar neutrinos that originate from 8 B decay. In the same figures we compared the 
results of our code BC (solid lines) with the approximate analytic results of Ram (dot- 
ted lines). (There are no dotted lines in the LR plots because the function f" B (E, e) 
has not been previously calculated.) As we mentioned earlier, Ram's formulae were 
computed in the small e approximation, keeping the logarithmically divergent terms 
proportional to ln(e/m), but neglecting the remaining e-dependent terms. Our results 
confirm Ram's ones in the small e limit. If e is not small, the discrepancy between 
solid and dotted curves increases with increasing values of e. Figures 1 and 2 also 
show that the dotted curves are not always positive. This is of course an unphysical 
property because the HB differential cross section, being a transition probability for a 
physical process, cannot be negative. (We also note that, in full generality, we could 
set g L or g R to zero, in which case the HB differential cross section would consist only 
of its R or L parts, respectively.) Our functions f BB (E,e) and f BB (E,e) are always 
positive (or zero). 

The total 0(a) QED corrections f L (E) and f R (E), given by the sum of V, SB 
and HB contributions (see eq. (0)), can be easily obtained by adding the analytic 
results of eqs. (@) and (H) to Ram's HB (lengthy) ones. Both SB and HB corrections 
were computed in the small e approximation, and the logarithmically divergent terms 
proportional to ln(e/m) exactly drop out upon adding these soft and hard contribu- 
tions. The remaining e-dependent terms, which were neglected in both SB and HB 
calculations, must cancel in the sum as well, and are therefore irrelevant in the com- 
putation of the total QED corrections of eq. ([]). The LR case is slightly different: 
Ram's formulae, which were used to derive the small e approximation for f" B (E,e) 
and f R B (E,e), do not provide us with the corresponding LR correction. In order to 
compute f LR (E) we have therefore added the V and SB analytic results of eqs. (|2|) 
and @ to our HB numerical results. The "exact" e dependence of our HB results is 
not completely canceled by that of the SB, which includes only terms proportional 
to ln(e/m), and the sum f LR (E) contains therefore a residual (not logarithmically 
divergent) dependence on the photon energy threshold e. This spurious dependence 
has been minimized by fixing e to be a very small value e LR chosen so as to have an 
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estimated induced relative error as small as 0(0.1%). Q 

Equations (|2|) and (|6|) determine the analytic expression of f^ s (E, e) (the infrared- 
finite sum of V and SB corrections) in the small e approximation. But the complete e 
dependence of our numerical f BB (E, e) computations, combined with the knowledge of 
the above described f x (E) functions, allows us to determine also the "exact" f^ s (E, e) 
functions via the subtraction 



ft s {E,e) = f x {E)-f BB {E,e). (12) 

These will be the "exact" VS corrections employed in the rest of our analysis. 

In fig. 3 we plotted the functions f x {E) (thick solid), f" B (E,e) (medium solid) 
and f% s (E,e) (thin solid) for v = 0.862 MeV. The threshold e in the VS and HB 
functions was set to 0.02 MeV. In figs. 4 and 5 we plotted the same functions with 
v = 10 MeV (e = 1 MeV) and v = 1 GeV (e = 50 MeV), respectively. 

In figs. 3, 4 and 5 we also plotted the simple approximate formulae for f x (E) in- 
troduced in ref. [16|] (dotted lines). These compact analytic expressions were obtained 
by modifying the expressions of ref. JTIJ, which had been evaluated in the extreme 
relativistic approximation. (The LR term of the differential cross section, being pro- 
portional to (m/z/), vanishes in the extreme relativistic limit and, therefore, cannot be 



derived from ref. ||10|| . As a consequence, the LR approximation of ref. |16| is only a 
(very educated!) guess.) Thanks to their simplicity, the compact formulae of ref. [|T(J 
are easy to use and are employed, for example, by the Super-Kamiokande collabora- 
tion [^T|] in their Monte Carlo simulations for the analysis of the solar neutrino energy 
spectrum. 



As it was noted in refs. @,[16||, a U fx(E) functions contain a term which diverges 
logarithmically at the end of the spectrum. This feature, related to the infrared di- 
vergence, is similar to the one encountered in the QED corrections to the /z-decay 
spectrum |f22],|23|]. If E gets very close to the endpoint we have {a/n)f x {E) « — 1, 
clearly indicating a breakdown of the perturbative expansion and the need to con- 
sider multiple-photon emission. However, this divergence can be easily removed (in 
agreement with the KLN theorem [[23|,[24]]) by integrating the differential cross section 
over small energy intervals corresponding to the experimental energy resolution. We 
also note that the singularity of f LR (E) foiE = m does not pose a problem, as the 
product [mz j v) f LR {E\ which appears in the LR part of the differential cross section, 
is finite in the limit E — > m. This can be seen from the dashed line in the LR plot of 
fig. 3, which indicates the function (mz/v)f hR {E). In the same plot, the dot-dashed 
line is the product of the f LR (E) approximation of ref. [pll and (mz/u). 



3 With the exception of E belonging to a tiny interval of 0(e LR ) at the endpoint -E m ax- We remind 
the reader that 0.1% is also the relative numerical uncertainty used by our code BC in the computation 
of the functions f" B (E, e) and produces a totally negligible relative error {a / if) f BB (E ^ e) x 0.1% in 
the corresponding X parts of the differential cross section in eq. (|]). 
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4 Spectrum of the Combined Energy of Electron 
and Photon 



We will now turn our attention to the analysis of the differential v\ + e — > v\ + e (+7) 
cross section relevant to experiments measuring the total combined energy of the recoil 
electron and a possible accompanying photon emitted in the scattering process. We 
will begin by considering bremsstrahlung events with a photon of energy uj larger than 
the usual threshold e (HB). 

The HB differential cross section [da/d(E + u;)]hb can be immediately derived 
from the HB corrections to the energy spectrum of the final neutrino. In the elastic 
reaction vi+e — > vi+e, the final neutrino energy v' ranges from z/' nin = vmj (2u+m) to 
^max = v (^ ne value v' = u' mhl occurs when the final electron and neutrino are scattered 
back to back, with the electron moving in the forward direction with E = Em&x, the 



value v' 



v' 



occurs in the forward scattering situation). When a photon of energy 



uj > e is emitted, v' varies between and v — e. If we now define the HB functions 



h(v',e, u) 



da 



and h(E + u,e,u) 
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da 
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such that 



h(i/, e, v) dv' 



v+m 



10 Jm+e 

conservation of energy implies then 



h(E + u,e, v) d(E + 



h(E + uj, e, v) = h(v' , e, u) = h{y + m — (E + cu), e, u), 

and the HB differential cross section with respect to the sum of the electron and 
photon energies can be directly obtained from [da/du']^- The variable E + uj varies 
between m + e and m + v (note that m + v = E max + f„ lin ). 

The function [da/dv']^, computed by our code BC, has been evaluated in a manner 
similar to the HB corrections to the electron recoil energy spectrum (see sect. 3). For a 
given initial neutrino energy, the five independent variables of the three-body phase- 
space describing the final state have been chosen to be the four angular variables of 
the final neutrino and photon, plus the final neutrino energy. In analogy with the case 
of the electron spectrum, we imposed the condition uj > e (once again, the threshold 
e is not assumed to be small and can vary up to the kinematic limit v) and evaluated 
the corresponding bounds on the chosen kinematic variables. Just as for the f" B (E, e) 
functions of sect. 3, the last integrations over the angular variables were performed 
numerically using the Monte Carlo method and requiring a very high (0.1%) relative 
accuracy. 

A check of the consistency of our results was performed by comparing the values of 
the total HB cross section ohb^, e) obtained by integrating the differential HB cross 



8 



sections of sects. 3 and 4. The equality 




dE 



da 
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d(E + u) 
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J HB 
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J HB 



has been tested for several values of v and e, and all relative deviations were found to 
be smaller than 0.1% (which is also the relative accuracy of the integrands). 

We now combine virtual, soft and hard bremsstrahlung contributions in order to 
evaluate the complete 0(a) QED prediction for the differential cross section da/d(E+ 
u) of reaction (||). In sect. 3 we computed the total QED corrections by simply adding 
these three parts. Their sum does not depend on the threshold e. The combination of 
VS and HB terms requires here a more careful analysis. We proceeded as follows. Let's 
consider an experimental setup for v—e scattering able to measure the photon energy 
if it's higher than a threshold e, but completely blind to low energy photons (uj<e). 
Let's also assume that the electron energy E is precisely measurable independently 
of its value. This detector can therefore measure the usual electron spectrum da/dE 
as well as the differential cross section da/dE^, where the variable E u is defined as 
follows, 



Figure 6 shows the E-u plane for E £ [m, m + u] and u £ [0, v\ (E, however, cannot 
exceed its elastic endpoint E max ). For a moment we will assume (as we did in drawing 
fig. 6) that e is smaller than e* = u' min , but we will later on free our analysis from 
this constraint. The vertical segments OP and PQ indicate, respectively, the sets 
of points contributing to the VS and HB corrections to the electron spectrum (see 
sect. 3) for a specific value E = E (in the virtual corrections it's simply u = 0). 
The overall QED corrections to the electron spectrum are obtained by adding up the 
VS and HB terms and clearly do not depend on the value of e. In fig. 6, diagonal 
lines indicate sets of points with the same value of the combined electron-photon 
energy. The set of all points in the E-u plane having the same E u value consists 
of one or two segments, according to the magnitude of E u which varies between m 
and m + u. If m < E^ < m + e, the photon energy is lower than the threshold 
e and the measured combined energy E u equals E. The QED corrections to the 
differential cross section da/dE^ will coincide, in this case, with the VS corrections to 
the electron spectrum da/dE (see e.g. segment AB). If m + e < E w < E max , the QED 
prediction for da/dE^ will consist of two contributions, namely the HB cross section 
[da/d(E + cj)]hb (see e.g. segment CD) plus the VS corrections to da/dE (segment 
EF). Finally, if E max < E^ < m + u, da/dE^ will be identical to [da /d ( E + u)]ub (see 
e.g. segment GH). It is important to notice that the complete QED prediction for the 
differential cross section with respect to E u depends on the threshold e, contrary to 
the complete QED prediction for the electron spectrum computed in sect. 3. 




(13) 



9 



We can summarize our results in a very simple way. The SM prediction for the 
spectrum of the combined energy of electron and photon in reaction (|8]) can be cast, 
up to corrections of 0(a), in the form 
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(14) 



where z u = (E u — m)jv and 9 = 9(E max — Ej). As we mentioned in sect. 3, the 
deviations of the functions g L (Eu) and g^E^) from the lowest-order values g L and 
g R reflect the effect of the electroweak corrections (for virtual corrections it is u — 
and E u — E). The functions f x (E i0 ,e,u) (X = L,R or LR), defined in the range 
[m,m + u], describe the QED effects discussed earlier in this section (once again, 
for simplicity of notation, we will drop their v dependence). Following our previous 
analysis, these functions can be written in the very simple form 



(15) 



where f^ 8 (E u , e) are the "exact" VS corrections of sect. 3 (eq. fll2|)) and the functions 
f" B (E^,e) are derived by dividing the L, R and LR parts of the above-studied HB 
cross section [da / d(E + uj)] hb by Cgl, Cg^l — zJ) 2 and —Cg L g R (mZu/v) respectively, 
with C = ImG^a I ix 1 . The 9 functions in eq. (0) reflect the fact that the lowest order 
prediction for da/ dE u has a step at E u = E max and is zero if E u lies outside the elastic 
range [m, E m3X ]. The VS functions /^(E^, e) are proportional to the same 9 function, 
while the corrections /^(E^, e) are set to zero if E w £ [m + e,m + v\. 

Earlier in this section we assumed e < e* (e* = v' min ). Nevertheless, the < e < v 
case can be discussed analogously and our simple prescription ( |15|) is valid in both 
cases. 

In fig. 7 we show, as an example, the function f L (E u , e) (thick) and its components 
/f s (^,e) (medium) and fl s (E^,e) (thin) for v = 0.862 MeV. The solid (dashed) 
lines correspond to the choice e = 0.2 MeV (e = 0.05 MeV). Once again we would 
like to emphasize the e dependence of the complete QED corrections fxiE^, e), to be 
contrasted with the e independence of the f x (E) functions of sect. 3. 

In figs. 8, 9 and 10 we compare the results of sects. 3 and 4. In fig. 8 we chose 
v = 0.862 MeV and plotted the functions f x (E) (thick) and J x (E m e) for e = 0.1 MeV 
(medium) and 0.001 MeV (thin). In figs. 9 and 10 we plotted the same functions with 
v = 10 MeV (e = 1 MeV, 0.1 MeV) and v = 1 GeV (e = 100 MeV, 10 MeV), 
respectively. 

We would like to remind the reader that the functions f x (E) can be obtained from 
f x (E UJ , e) by simply setting e = v. The limiting case e = was studied in detail in 
ref. |12] (in particular, the results of the second article of this reference were obtained, 
like ours, without employing the ultrarelativistic approximation E ^> m). 
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5 Discussion and Conclusions 



When are the results of sects. 2, 3 and 4 applicable? In sects. 2 and 3 we evaluated 
the 0(a) SM prediction for the electron spectrum in the reaction v\ + e — ► P\ + e (+7) 
(eq. (|J)), where (+7) indicates the possible emission of a photon. In this calculation 
we assumed that the final-state photon is not detected and, as a consequence, we 
integrated over all possible values of the photon energy uj. Therefore, eq. (|9]) is the 
appropriate theoretical prediction to use in the analysis of v—e scattering when the 
detector is completely blind to photons of all energies, but can precisely measure E, 
the energy of the electron. Of course, a detector could provide more information by 
detecting photons as soon as their energy is above an experimental threshold e. In 
this case, still assuming a precise determination of E, one can employ eq. (j^), minus 
its HB correction, to analyze those events which are counted as nonradiative (elastic), 
while the HB part can be used, at least in principle, for a separate determination of 
the inelastic cross section. Indeed, contrary to previous calculations, our predictions 
are valid for an arbitrary value of the threshold e (and include the previously unknown 
LR term). 

In sect. 4 we evaluated the spectrum of the total combined energy of the recoil elec- 
tron and a possible accompanying photon emitted in the scattering process (eq. fll4]); 
E w was defined in eq. ([13])). This type of analysis is useful when the photon energy 
uj cannot be separately determined although it fully contributes to the precise total 
energy measurement if its value is above a specific threshold e. Let's consider an exper- 
imental setup able to measure the photon energy if it's higher than e, but completely 
blind to low energy photons (u<e). Let's also assume that the electron energy E is 
precisely measurable independently of its value. This detector can determine both the 
differential cross section dcr/dE^ (eq. ([14])) and the electron spectrum da/dE (eq. ©) 
(as well as its separate HB component). There are experiments, however, which can- 
not measure E, but only E w , with a specific value of the threshold e. BOREXINO [p5fl 
and KamLAND pBf , for example, are liquid scintillation detectors in which photons 
and electrons induce practically the same response. If a photon is emitted in the v—e 
scattering process, its energy uj is counted together with E, provided their sum lies 
within a specific range. The appropriate theoretical prediction for their analysis is 
given, therefore, by the cross section da/dE w of eq. ([L4]) with a very small value of e 
(for the case v = 0.862 MeV see the thin lines in fig. 8). However, we should point 
out that although the QED corrections f x (E) (in eq. (|9])) and f x (E UJ ,e) with small e 
(in eq. (JT4])) are different, their numerical values are very small when v = 0.862 MeV, 
the energy of the monochromatic neutrinos produced by electron capture on 7 Be in 
the solar interior. In fact, as shown in fig. 8, both (a/ii)f x (E) and (a/^fx^E^, e) 
with small e are in this case of 0(<1%), and neither of the above collaborations is 
likely to reach this high level of accuracy in their analyses of the crucial 7 Be line. 

There are detectors in which it might not be possible to identify the measured 
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energy with either E or E u . Indeed, the electron and the photon may produce indis- 
tinguishable signals and the total observed energy might not be the simple sum of E 
and uj, but some other function of these two variables. Super-Kamiokande (SK) |21|| , 
for example, a water Cherenkov counter measuring the light emitted by electrons 
recoiling from neutrino scattering, uses the number of hit photomultiplier tubes to 
determine the electron energy. However, a photon emitted in the scattering process 
may induce additional hits indistinguishable from those of the electron. Moreover, 
a photon and an electron of the same energy may produce different numbers of hits 
and, therefore, it might not be possible to identify the total measured energy with the 
sum E + uj. 

SK measures solar neutrinos with energies varying from 5 to 18 MeV. For v = 
10 MeV, fig. 9 shows that the QED corrections to the differential cross sections da/dE 
(eq. (H)) and da /dE u (eq. (|T4])) are of 0(1%). Corrections of this order may be relevant 
for the analysis of the very precise data obtained by this collaboration. In fact, SK's 
Monte Carlo simulations of the expected energy spectrum of recoil electrons from solar 
neutrino scattering include the QED corrections of ref. |16|] (as well as the EW ones). 
As we investigated in sect. 3, these corrections provide good approximations of the 
complete 0(a) QED corrections f x (E) to the electron spectrum of eq. (||) (see fig. 4). 
Our previous discussion, however, seems to suggest that these corrections may not be 
appropriate for SK's solar neutrino analysis. On the other hand, the SM prediction for 
the spectrum of the combined energy of electron and photon of sect. 4 (eq. (0)) may 
be suitable, but only if we can assume a similar efficiency in the detection of photons 
and electrons, and if also relatively low energy electrons contribute to the total energy 
measurement. If these conditions are not met, and the precision of the data requires 
it, one should probably perform a dedicated analysis of the double differential cross 
section d 2 a / \dE duj) with a response function specifically designed for this detector. 
A triple differential cross section d 3 a/(dE du d(p), where is the angle between the 
directions of the electron and the photon, may also be useful (see the third article of 
ref. [|, and ref. 0). 

We will conclude by noticing that the QED corrections, contrary to the EW ones, 
depend strongly on the initial neutrino energy and become sizeable at high energies. 
Their appropriate expression will have to be taken into account in the analysis of 
future precise high energy u—e scattering experiments. 
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Figure 1: Comparison of the f" B (E, e) functions (solid lines) with their small e analytic 
approximations (dotted lines) for an initial neutrino energy v = 0.862 MeV. Thick, normal 
and thin lines indicate, respectively, e = 0.001, 0.01 and 0.1 MeV. 
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Figure 2: Comparison of the f" B (E, e) functions (solid lines) with their small e analytic 
approximations (dotted lines) for an initial neutrino energy v = 10 MeV. Thick, normal 
and thin lines indicate, respectively, e = 0.5, 1 and 2 MeV. 
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Figure 3: The functions f x (E) (thick solid), f" B (E, e) (medium solid) and f^ s (E, e) 
(thin solid) for v = 0.862 MeV and e = 0.02 MeV. The dotted lines represent the 
f x (E) approximations of ref. flTjJ. In the LR figure, the dot-dashed line is the product 
of the f LR (E) approximation of ref. |TB[ and (mz/v), while the dashed line indicates 
the product {mz/u)f LR {E). 
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Figure 4: Same as Fig. 3, but for v = 10 MeV and e = 1 MeV. The dashed and 
dot-dashed lines are very close to zero and are not indicated. 
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Figure 5: Same as Fig. 3, but for v =1 GeV and e = 50 MeV. The dashed and 
dot-dashed lines are very close to zero and are not indicated. 
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Figure 6: The E-u> plane for E G [m, m + u] and u> G [0, u\. See text for details. 
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Figure 7: The function f L (E^,e) (thick) and its components f" B (E LJ ,e) (medium) 
and f^ s (E UJ ,e) (thin) for v = 0.862 MeV. Solid and dashed lines indicate, respectively, 
e = 0.2 and 0.05 MeV. 
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Figure 8: The functions f x (E) (thick) and f x {E UJ ,e) for e = 0.1 MeV (medium) and 
0.001 MeV (thin), v = 0.862 MeV. 
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Figure 9: The functions f x (E) (thick) and f x (E UJ ,e) for e = 1 MeV (medium) and 
0.1 MeV (thin), v = 10 MeV. 
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Figure 10: The functions f x (E) (thick) and f x (E UJ ,e) for e = 100 MeV (medium) 
and 10 MeV (thin), v = 1 GeV. 
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